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Renormalization constants for multiplicatively renormalizable parity-odd four-fermion operators are computed in 
various different Schrodinger Functional (SF) schemes and lattice regularizations with Wilson quarks at one-loop 
order in perturbation theory. Our results are used in the calculation of their NLO anomalous dimensions, through 
matching to continuum schemes. They also enable a comparison of the two-loop perturbative RG running to the 
previously obtained nonperturbative one in the region of small renormalized coupling. 



1. Introduction 

Four-fermion operators occur in the transition 
amplitudes of many physical processes within the 
SM and beyond. An example is provided by the 
oscillations of the K'^ — system, where the 
mixing amplitude is given by a matrix element of 
the form 



(1) 



Renormalizing the operator Oyy^^^ on the lat- 
tice with Wilson quarks requires special care, as 
it mixes with other four-fermion operators. The 
situation is simpler in the framework of tmQCD, 
where the operator in question can be mapped 
onto its parity-odd counterpart OyA-i-AV' which 
is multiplicatively renormalizable "7 . Nonpertur- 
bative renormalization of the latter in SF schemes 
has been reported in where its scale depen- 
dence has been computed in the continuum limit 
through finite size scaling techniques. On the 
other hand, the corresponding one-loop pertur- 
bative calculation is required in order to make 
contact to continuum schemes. It can also be 
used for an analytical estimation of the perturba- 
tive lattice artefacts. 



2. SF renormalization 

Renormalization of a local fermion operator is 
pursued in a SF scheme by imposing a suitable 
renormalization condition at vanishing renormal- 
ized quark mass on a correlation function, prop- 
erly defined to probe the operator by bilinear 
boundary quark sources 



%/2[r] 



'E C/,(y)rC/.(z) 



s'M.n = C},(y')rc;.(z'). (2) 

y',z' 

Here the boundary fields C, and C, are defined at 
xq = 0, while C' and C' are defined at xq = T; F is 
a generic Dirac matrix and fi are flavour indices. 
The simplest way to probe the four-fermion op- 
erator 



^VA+AV 



= \ [V'17mV'2V'37m75'04 + 

V'l7Ai75V'2V'37MV'4 ± (2 ^ 4)] (3) 



*Talk given at Lattice 2004 by F.P. 



is by inserting it in a correlation function with 
three bilinear boundary sources 

'^^(^o) = ^{S'^A^:i]0^{x)S2i[Ti]S,,[T2]). (4) 

Parity conservation implies that one has to pro- 
vide a parity-odd choice of the set [Fi,F2,r3] in 



1 



2 



order to get a nonvanishing correlation function. 
Provided this, the uhraviolet divergences of the 
boundary fields must be removed in order to iso- 
late the divergences of the operator. To this pur- 
pose /i*(xo) can be normalized as 



/r t 



(5) 



Q+/3=3/2 

where we have introduced the standard boundary 
correlation functions with no operator insertion 

/, ^ ' 

fci = 



1 

3Z6 



(6) 



fe=i 



The condition a + (3 = 3/2 makes h^{xQ) free of 
boundary divergences and suitable to be renor- 
malized through a multiplicative condition 

h%{T/2) = ht^^{Tl2) atm^^O, (7) 

where the subtraction point is at /i = l/L, niji 
is the renormalized quark mass, the operator 
CvA-i-AV placed in the middle of the time extent 
of the SF and 

hi{Tl2) = Z±,+Av(3o, a/i) h^{Tl2). (8) 

Apart from the indicated dependence, the renor- 
malization constant ^va+av depends upon every 
detail, such as the value of T/L (in this work 
T/L = 1), the choice of the Dirac structure of 
the sources and the exponents a, (3 in eq. (5). A 
number of possible combinations exist, each one 
generating a different renormalization scheme in 
the framework of the SF. In Table 1 we report the 
ones we have considered. 

3. One loop perturbative expansion 

In order to extract the one-loop contribution to 
the renormalization constant, all the correlation 
functions must be expanded in powers of the bare 
coupling. 
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dmo 



0{gt), (9) 



with X being h^, /i, fci, as well as the renormal- 
ization constant itself. 



■^VA+AV ~ 1 



2;,±(1) 

i/O^VA+AV 



(10) 



scheme 



[ri,r2,r3] 



a 



f3 



I 


[75,75,75] 


3/2 





II 


6 Sijfc=l ^ijk [li , 7i , 7fe] 


3/2 





III 


iELi[75,7fc,7fc] 
5ELi[7fc,75,7fe] 
iELi[7fc,7fc,75] 


3/2 





IV 


3/2 





V 


3/2 





VI 







3/2 


VII 


5ELi[75,7fc,7fe] 
5ELi[7fc,75,7fe] 
iELi[7fc,7fc,75] 


1/2 


1 


VIII 


1/2 


1 


IX 


1/2 


1 



Table 1 

Terms proportional to the one-loop coefficient 
of the critical mass mi^'' are needed to set the 
correlation functions to zero renormalized quark 
mass. The values of f[^'^^ and k^'^"^ are known 
from the literature [S], while have been 

computed through numerical integration of the 
related Feynman diagrams for both Wilson and 
Clover actions. 

The analytic structure of ^yl+AV provided by 



the "log" -divergent relation 

7±(1) _ „± , ^+(0)1^ 

^VA+AV - ^SF +7 ' ' in - 



L 



O 



where = ±6(1 T l/N)/{4:n)^ is the LO 

anomalous dimension of CyA+AV ^^d B^p is a 
finite term, which depends upon the renormaliza- 
tion scheme and can be computed by inserting the 
perturbative expansion of the correlation func- 
tions and the renormalization constant in eq. (7). 



4. Matching to DRED 

Given the NLO (scheme dependent) anomalous 
dimension in a reference scheme, it can be com- 
puted in any other scheme, provided the one-loop 
order renormalization constant is known in the 
latter. In our case we choose DRED as the refer- 
ence scheme, and the matching formula is given 

by 



±(1) 

7SF 



±(1) 

7dred 



2bo 



+ 7^ 



±(0), 



(12) 
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-7+ 

1 


(1) X 10''' 


I 


0.517(3) 


+ 0.2032(2) Nt 


II 


-6.072(5; 


+ 0.6026(3) Nf 


III 


1.690(3) 


+ 0.1321(2) Nf 


IV 


-7.966(5' 


+ 0.7174(3) Nf 


V 


-7.290(5' 


+ 0.6763(3) Nf 


VI 


-7.830(8] 


+ 0.7090(5) Nf 


VII 


0.519(5) 


+ 0.2031(3) Nf 


VIII 


-9.139(8; 


+ 0.7885(5) Nf 


IX 


-8.463(8; 


+ 0.7475(5) Nf 


scheme 


7" 


(1) X 102 


I 


2.3660(3) 


- 0.19711(2) Nf 


II 


0.9320(5) 


- 0.11020(3) Nf 


III 


2.4833(3) 


- 0.20422(2) Nf 


IV 


0.5530(4) 


- 0.08723(3) Nf 


V 


0.8057(4) 


- 0.10255(3) Nf 


VI 


0.7559(5) 


- 0.09953(3) Nf 


VII 


2.3660(3) 


- 0.19711(2) Nf 


VIII 


0.4355(4) 


- 0.08011(3) Nf 


IX 


0.6882(4) 


- 0.09543(3) Nf 



Table 2 



where the second term on the RHS measures the 
shift of the finite term of the one-loop renormal- 
ization constant and the third one measures the 
shift of the coupling between the schemes. The 
and Cg can be retrieved 



coefficients 7d^eD' 



from the literature |4I5I6I7| and 6o is the LO con- 
tribution to the /3-function. In Table 2 we report 
the (preliminary) values of 7='=(^) in the schemes 
defined in Table 1. The dependence of the NLO- 
AD upon Nf is implicit in the coefficients 7dreD' 
6o and Cf. 



5. NP vs. NLO running of O 



± 

VA+AV 



The knowledge of the NLO anomalous dimension 
allows, in particular, to compare the nonpertur- 
bative and NLO (scheme dependent) running. To 
do this, we consider the step scaling function (ssf) 



' VA+AV V 



lim 



Z. 



± 

VA+AV 



(5o,a/2i) 



*o 2'vA+Av(5o,a/i) 



(13) 



which has been computed nonperturbatively in 
and expand it to NLO: 

^VA+Av(") = 1 + + ^2^ + 0{u^) 



at - 7**°^ In 2, 
at = 7^<^^ 

rl 



In 2 



2(7^(°))^ + 6o7^(°)](ln2)= 



(14) 



As an example, in Figure 1 we report plots of 



the ssf a 



VA+AV 



in the schemes I and II. In both 



plots, the solid line represents the LO running, 
the dashed line represents the NLO running and 
points are the nonperturbative data. The agree- 
ment between perturbation theory and nonper- 
turbative data is scheme dependent; in some 
schemes the agreement improves when switching 
to the NLO, while in some others it worsens. 




Figure 1. o"va+av(*') (schemes I and II): LO running 
(continuous line) , NLO running (dashed line) and NP 
data (points) from pp. 

We conclude that the SF scheme has to be chosen 
with care in order to achieve a good control of the 
matching to the continuum schemes. 
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